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We obtain exact analytic expressions for a class of functions expressed as integrals over the Haar 
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from ensemble theory. We further generalize these results to arbitrary eigenvalue distributions, 
allowing a detailed comparison of typical regular and chaotic systems with the help of concepts from 
random matrix theory. To illustrate the physical relevance and the general applicability of our results 
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dynamical maps, the generic dynamics of system-environment pure state entanglement and, finally, 
the equilibration of generic open and closed quantum systems. 
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I. INTRODUCTION 

The theory of open quantum systems allows for the de- 
scription of quantum systems interacting with their en- 
vironment [1]. Recently, the focus of research has been 
drawn towards more and more complex models and even 
biological systems have been approached by the methods 
of open systems [2-4]. Increasing size and complexity 
lead to more interesting and realistic models but often 
cause the treatment of the problem to become very cum- 
bersome, if not impossible. In this paper, instead of in- 
vestigating a special system, we are looking for generic 
features which emerge in a variety of complex systems. 
We study these features employing a statistical approach 
based on averages over a multitude of realizations as it is 
usual practice in ensemble theory for statistical physics 
[5]. Similar approaches have already been applied to 
complex quantum systems [6-12]. We analyze dynamical 
features of open systems by averaging over the unitary 
time evolution operator, which determines the dynamics 
of system and environment. The thereby obtained results 
reflect the behavior of a generic complex open quantum 
system. 

The average values are expressed as integrals employ- 
ing the uniform Haar measure on the unitary group. All 
integrals of interest in this paper can be written as func- 
tions of the form 

EW(Xi,...,X n _i) - j dv(U)tfX 1 U...tfX n - 1 U, 

(1) 

where X\,...,X n —i are arbitrary operators on a d- 
dimensional Hilbert space T-L and d/j,(U) denotes the Haar 
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measure on the unitary group U{d). The expression on 
the right-hand side of Eq. (1) involves an integration 
over n unitary matrices and in the following the func- 
tion E^(Xi, . . . , A"„_i) will be called the nth moment 
function of the unitary group. The results of Collins and 
Sniady [13] allow, in principle, for an exact determination 
of these moments up to arbitrary order. In this paper, 
moments up to the order of n — 8 will appear. For the 
sake of readability, the explicit derivation of these expres- 
sions is shifted to the appendix. 

Equipped with these mathematical tools we are able 
to analytically tackle a variety of physically motivated 
expressions. In Sec. II, we present three general results 
which will form the foundation for most of the following 
analysis in this paper. These results involve the expec- 
tation value of the norm of an arbitrary operator, which 
has been subjected to a unitary time evolution before the 
environment is traced out. We obtain a uniform average, 
as well as the corresponding variance. Additionally, a 
decomposition of the Hamiltonian into eigenvalues and 
eigenvectors allows to distinguish between the generic ef- 
fects of typical systems with well-known energy-level dis- 
tributions, which will be put into context with random 
matrix theory and quantum chaos [9, 10, 14, 15]. 

In subsequent sections, the versatility of the general 
results is demonstrated by a study of various applica- 
tions which are of current interest in the fields of open 
quantum systems [1], quantum information theory [16], 
quantum chaos [14, 15] and quantum thermodynamics 
[17]. In Sec. Ill we investigate the generic signature 
of system-environment correlations in the open-system 
dynamics. In Sec. IV we identify the average reduced 
dynamics of an open quantum system as the depolariz- 
ing quantum channel. Section V deals with the average 
system-environment pure state entanglement dynamics. 
Finally, in Sec. VI, we apply our approach to the ther- 
malization process of generic complex quantum systems, 
examining both isolated and open systems. 



2 



II. AVERAGE HILBERT-SCHMIDT NORM 
EVOLUTION IN AN OPEN QUANTUM SYSTEM 

We consider an open quantum system, composed of 
the Hilbert spaces of system Hs and environment He 
with dimensions d$ and d Bl respectively. The com- 
position of system and environment is assumed to be 
closed, such that we obtain unitary dynamics for the to- 
tal system. Furthermore, we denote by M an arbitrary, 
fixed self-adjoint operator with spectral decomposition 
M = yi —1 rrii \i){i\, acting on the d-dimensional Hilbert 
space H = Hs <8> He of the closed system. In the fol- 
lowing we are interested in the open system part of M 
after being exposed to the unitary time evolution of the 
total system. More precisely, we introduce the following 
operator: 



A* = Ti B {U t MU}}, 



(2) 



where Ut denotes the unitary time evolution operator, 
propagating the total system states from time to t, and 
Tte the partial trace over the environment. 

Our goal in this section is to obtain an expectation 
value of the squared Hilbert-Schmidt norm of A t , av- 
eraging over all unitary time evolution operators U t . 
We choose the squared Hilbert-Schmidt norm ||A|| 2 = 
Tr{A^A} since it enables us to analytically obtain ex- 
act expressions for the expectation values in terms of the 
functions E'"' defined in Eq. (1). 



A. The uniform average 

In our first approach we replace the time evolution op- 
erator Ut by a random matrix U and average over all pos- 
sible realizations of unitary matrices employing the Haar 
measure d/j,(U). This is an invariant measure, giving uni- 
form statistical weights to all unitary matrices. We use 
the following notation for operator-valued averages over 
the Haar measure: 



(F(U)) = / dfi(U)F(U). 



(3) 



This approach leads to time-independent results as the 
time argument is lost in the averaging process. Later we 
generalize this approach in a way which allows to main- 
tain the explicit time dependence of the average value. 



1, The expectation value 

We now derive the unitary average of the function 
||A|| 2 = \\T, E {UMW}\\\ Let {\ Vi ))ti and {\x,)}% 
be fixed orthonormal bases oiHs and He, respectively. 



The elements of the matrix A can be expressed as 
A fe/ =^m, (ip k Xj\U\i)(i\U^ \<piXj) 

= ^m i (i\U i A k iU\i), (4) 

i 

where we have defined the operators Aki = \(pi)(<pk\ ® I 
with the property 4i = Aik, I denoting the identity on 
He- The squared Hilbert-Schmidt norm of the Hermitian 
matrix A then reads 

||A|| 2 = TrA 2 = ]TA fei A ifc 

= EE m * m i w tf A * u w oi uiA lu ii) ■ (5) 

The average value of this expression can be written in 
terms of the general fourth moment function E^ 4 ) of the 
unitary group, given in the appendix in Eq. (A8), with 



X 1 - X- 



A k i and X 2 



(iiaii 2 ) = EE m ^ m {4) (A kl , ,4,) ii) . (6) 

k,l i,j 

Using d = dsdE and the relations Tr{4 ;J 4fc;} — &e and 
TrA k i = Tr4 ; = 5 k id E , we obtain the general result for 
the uniform average [10]: 



Tr £ {UMU^}\\ 2 ) = Ci \\M\\ 2 + C 2 (TrAf) 2 , (7) 



with 



Ci 



d 2 s d E -d E , „ d s d 2 E - d s 
and Co ~ 



d\d\ 1 



44 - 1 



(8) 



Note that in the limit of high environmental Hilbert 
space dimensions, the coefficient C\ vanishes while C 2 
approaches the asymptotic value of 



2. The variance 

Next we derive the corresponding variance of the av- 
erage value (7), which is defined as 



Var(||A|| 2 ) = (||A|| 4 )-<||A|| 2 ) 5 



(9) 



Since the average value is already known, the task is to 
determine the quantity (||A|| 4 ). In terms of the matrix 
elements given in Eq. (4) this yields 



|A|| 4 - E AyA^AjyA 

i,j,k,l 



Ik 



771 a HI ^llljnis (a\ AijU | 
i : j,k,l 

x (/3|[/U Jl (7|/3)(7|C/U fci C/| 7 )((5|C/ t J 4 ifc C/|(S). 

(10) 
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Here, eight unitaries are involved in the integration. 
Hence, we can obtain the result with the aid of the gen- 
eral eighth moment function E^ 8 ) , given in the appendix 
in Eq. (All), making the following choice of matrices: 

Xi = X\ = Ajj, X 5 = X 7 = Am, 

X 2 = \a)(fi\, X 4 = |/?>< 7 |, X e = \ 7 ){6\. (11) 

The variance is given according to Eq. (9), 

Var(|| A|| 2 ) = ci(TrM) 4 + c 2 (TrM)(TrM 3 ) 

+ c 3 (TrM) 2 (TrAf 2 ) + c 4 (TrM 2 ) 2 

+ c 5 (TrM 4 ), (12) 

with the coefficients 

ci = 26(11 + 44), 
c 2 = 40a, 

c 3 = -4bd s d E {ll + d 2 s d%), 

c 4 = 26(15 - 444 + d 4 s d%), 

c 5 = -l0ad s d E , (13) 

where 

a = (4 - 1)(4 - 1)(44(44 - 7) 2 - 36)"\ 

and 

b = (4-i)(4-i)(44-i)~ 2 (36-i344+44)~ 1 . 

For increasing environmental dimension d E , all coeffi- 
cients ci, . . . , C5 tend towards zero, indicating a general 
universality of the average value (7) for high Hilbcrt 
space dimensions, see Fig. 1. This result is in agreement 
with the phenomenon of the concentration of measure, 
which is based on Levy's lemma [18, 19]. The decreas- 
ing variance for this class of average values also endorses 
the statement that high-dimensional generic complex sys- 
tems are well described by the average value. 

B. The average for general eigenvalue distributions 

1. The general result 

As a generalization of the result of the previous sec- 
tion, we decompose the Hamiltonian into eigenvectors 
and eigenvalues as H — WDW\ where W denotes a 
unitary matrix containing the eigenvectors of H and 
D = diag(£'i, E2, ■ ■ ■ ) represents a diagonal matrix with 
the eigenvalues of H as diagonal elements [9]. The uni- 
tary time evolution operator can be written as U t — 
We~~ lDt W^ . Assuming the eigenvectors to be random 
and independent of the eigenvalue distribution, we inte- 
grate over W with the uniform Haar measure and thereby 
obtain the generalized version of Eq. (7): 

(i|A|| 2 ^ = j dn(W) ||Tr B {We~ iDt Wl ' MWe lDt W^}\\ 2 , 

(14) 




FIG. 1. (Color online) Coefficients of the variance for ds = 
2 as a function of d_E. All coefficients asymptotically tend 
towards zero in the limit of high environmental dimensions. 



where we have defined the system matrix A = 
Tt e {We- lDt W^MWe lDt W^}. The time-dependence 
has been retained within this approach in contrast to the 
ansatz of Sec. II A. The following procedure is however 
very similar to the previous derivation. In terms of the 
matrix elements of A, 

A kl = "V. (Ml We iDt W^ A kl We~ iDt W^ |/x) , (15) 

the average value (14) takes the form 
(||A|| 2 ) = ]T m^m v 

x ( M | E( 8 V Dt , A kl , e~ wt , \n){v\ , e iDt ,A lk , e~ im ) \v) . 

(16) 

We have made use of the general result of the eighth 
moment function E^ 8 \ see Eq. (All), with the special 
choice of the matrices 

X\ =X 3 = Xl =X 7 = e-* Dt , 

X 2 = xl=A kh X 4 = \ix){v\. (17) 

Finally, we obtain 

(||A|| 2 ) =c x (t) \\M\\ 2 + Z 2 {t) (TrM) 2 

+ c 3 (t) \\Tr E M \\ 2 +c 4 (t) \\Tt s M\\ 2 , (18) 

where the coefficients ci(f),. .. , c 4 (i) are functions of the 
dimensions of system and environment and the distribu- 
tion of the eigenvalues Ej of the Hamiltonian. Introduc- 
ing the Fourier transform of the level density 

/(*) = JiS^ 4 ' ( 19 ) 
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these coefficients can be expressed as: 



ck(t) = 

Mt) = 
Mt) = 

Ci(t) = 
where 



[(d 2 - 3d%)b(t) - 2d 2 (d 2 E - 3)n{f(t) 2 f*(2t)} 

+ (d 2 ~9)(d 2 -d%)]/(ad E ), 

[d(d 2 E - 3)6(t) - 2d(d 2 - 3d 2 E )dl{f(t) 2 r(2t)} 

+ d(d 2 -9)(d 2 E -l)]/(ad E ), 

- [(d 2 ~3)b(t)+4d 2 n{f(t) 2 r(2t)}]/a, 

[2db(t) + 2d{d 2 - 3)5ft{/(i) 2 .f (2t)}] /a, (20) 



a = tf 1 - 10d 2 + 9, 
6W=4|/(t)| 2 -|/(2i)| 2 -d 2 |/(i)| 4 , 

and 5ft denotes the real part of a complex number. We 
observe that the coefficients ci(t), . . . , 04(f) depend only 
on the dimensions of system and environment and on the 
eigenvalue distribution of the Hamiltonian via functions 
of/(t). 
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FIG. 2. (Color online) Comparison of the average coefficient 
(ci(i)) for regular (Poi) and chaotic (CUE) systems for ds = 2 
and different environmental dimensions. 



2. Averages for regular and chaotic systems 

From the theory of quantum chaos and random matrix 
theory it is well known that the energy distribution of 
complex systems determines whether the system dynam- 
ics exhibits regular or chaotic behaviour [14, 15]. In order 
to distinguish between different ensembles, we average 
the set of eigenvalues of the total system Hamiltonian ac- 
cording to well-known characteristic statistics. Since the 
average value (18) depends on the energy spectrum via 
the function /(f), the averages of |/(f)| 2 , $t{f(t) 2 f* (2f)}, 
and |/(f)| 4 with respect to the eigenvalue distributions in 
question are required. In Appendix B, we derive the av- 
erages of these functions for regular systems, being rep- 
resented by Poissonian distributed eigenvalue spacings 
(Poi), and for chaotic systems with no time reversal sym- 
metry as represented by the Gaussian unitary ensemble 
(GUE). 

As an example we show in Fig. 2 the average time 
evolution of the coefficient ci(f) for regular and chaotic 
systems. Note that here and in the following the time f is 
taken to be dimensionless, the relevant unit of time being 
of the order of the quantity H/AE, where AE represents 
the width of the level density R\{E) (see Appendix B). 



tend towards nonzero functions of time: 
0, 



lim (ci(t)) Po 

d,E—>oo 

lim (c 2 (f)) Po 

d,E— >oo 



128f 4 + 4cos(4£) - cos(8f) - 3 



128d s t 4 



lim (c 3 (t)) p i = 

dE—too 

lim (c 4 (t)) Poi 



cos(f) sin(t) 



0. 



(21) 



If we additionally take the limit for large times t, the 
coefficients C2 (t) and C3 (t) yield the following asymptotic 
values: 



lim lim (c 2 (t)) Poi = 

t^roo ds— too US 



lim lim (c 3 (f)) Po 

t— >oo (Ie— too 



0. 



(22) 



This coincides with the time-independent result of the 
uniform average in the same limit, see Eq. (7). 

Now we investigate the Gaussian unitary ensemble. 
For large d E and fixed ds we obtain the following results 
(see Appendix B 2): 



3. Asymptotic values of the coefficients 

The result of the average value (18) depends crucially 
on the dimensions of system and environment. In this 
section we analyze this dependence and present a con- 
nection to the uniform average, given in Eq. (7). 

For the Poissonian level spacing distribution the coef- 
ficients Ci(t) and Ci(t) disappear in the limit of large d E 
with a fixed value of ds while the other two coefficients 



lim (ci(t)) GU E 



= 0, 



lim (c 2 (t)) que 



lim (c 3 (i)) 



GUE 



lim (c 4 (£))gue 



,h{2t) 



= 0. 



Ji(2t) 



(23) 
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After the additional limit for large t, again, we find that 
the values coincide with the results which we obtained 
for the Poissonian average and the uniform average: 



lim lim (c 2 (t)> gue = -7- 

t— >oo d E — >oo as 

lim lim (c 3 (i)) que = 0. 

t— >CG dE—tCG 



(24) 



While for the Poissonian statistics the asymptotic values 
for t — > oo depend on ds and in a complicated way, 
here we find that within the approximation for large val- 
ues of d (see Appendix B 2), the Gaussian average is equal 
to the time-independent result of the uniform average: 



lim (ci(t)) G uE s 

t—*-oo 


d 2 s d E - 


- d,E 






' 44 


- 1 




lim (c 2 (*))gue s 

t— >oo 


d s d% - 

J 44 


-d s 
- 1 


= c 2 , 




lim (c 3 (t)) gue p 

t— >OQ 










lim (c 4 (i))GUE - 

t— too 


s 0. 






(25) 



C. The effective environmental dimension 

Although the dimension of the system Hilbert space is 
usually quite obvious, from a physical point of view, it 
is not always clear which is a suitable dimension for the 
environment. In many cases, the environment consists 
of a huge Hilbert space of which only a small part actu- 
ally affects the system dynamics. Consider a fixed sys- 
tem, composed of system and environment with Hilbert 
space H = Hs ® He- Formally, it is always possible to 
add another Hilbert space He' to the system by defin- 
ing a new environment He — He ffi He' ■ The new 
composite Hilbert space then reads H = Hs ®He = 
(Hs ®He) © {Hs ®He')- One could interpret this as 
the inclusion of a part of the environment in the calcu- 
lations which is by no means coupled to the system of 
interest. Clearly, the additional environment He' will 
not affect the system dynamics at all, although formally 
the environmental dimension has been increased. 

In conclusion, the effective environmental dimension, 
to be inserted as d# in the formula for the averages, is 
not always given by the mathematical dimension of the 
environmental Hilbert space, but rather characterizes a 
certain subspace of the Hilbert space which plays a role 
in the system's dynamics. For a given physical model the 
effective environmental dimension may be hard to obtain 
and further investigations are certainly necessary in or- 
der to propose a clean definition and an applicable and 
efficient way for its theoretical and experimental deter- 
mination [20]. 



III. CORRELATIONS IN BIPARTITE SYSTEMS 

An important special case of the previously derived 
general results is obtained if M = p — p', with two states 
p and p' of the closed system, is inserted into Eqs. (7) 
and (12). The uniform average yields [10] 



\Tr E {U(p-p')U^}\ 



d 2 s d E 



44 - 1 



Wp-p' 



(26) 



with the variance 



Var(||Tr B {U(p- f/)rf}\\ ) = c 4 (Tr{(p- p'f}f 

+ c b Tv{(p-p'Y}. (27) 

Distance measures are particularly interesting in the 
context of geometric measures of correlations [21]. The 
distance of a given state p to an associated uncorrelated 
state can be interpreted as a measure for the amount of 
correlations in p. A possible measure for the total corre- 
lations of a state p is the distance of the state p to the 
completely uncorrelated product state of its marginals 
p' = Ps®Pe, where ps = Tr^p and pe = Trsp denote the 
reduced density operators of p for system and environ- 
ment, respectively [22]. Choosing p' as the state which 
is produced from p by local dephasing in the eigenbasis 
of ps , one can quantify the amount of quantum discord 
contained in p [10, 20, 21, 23-25]. 

From Eq. (26), one can see the remarkable result that 
the average distance of the two reduced states ps(t) — 
Tv E {U t pUl} and p' s = Ti- E {U t p'Uj} is proportional to 
the distance of the initial states of the composite system. 
The above examples for p and p' then show that the total 
correlations or the quantum discord, given by the right- 
hand side of Eq. (26) have direct impact on the generic 
evolution in the open system. This can be used in order 
to detect properties of the initial states locally [10, 20]. 

This direct proportionality can also be observed for the 
average for general eigenvalue distributions whenever the 
marginal states of p and p' are equal. More precisely, 
from Eq. (18) we find 



|Tr B {We- iDt W\p - p')We lDt W^}\ 

= c!(£) Hp-p'II 2 , 



(28) 



if Tt E {p-p'} = and Tt s {p-p'} = 0. This is obviously 
the case if p' is chosen to be the corresponding product 
state of p, but also holds for the locally dephased refer- 
ence state [10, 20], as described before. The function C\ (t) 
is plotted in Fig. 2 for typical regular and chaotic sys- 
tems. Equations (21) and (23) show that in both cases, 
this function approaches zero for large d^. The same 
holds for the coefficient in the uniform case, displayed in 
Eq. (26). Hence, we find that the average evolution of an 
arbitrary pair of quantum states in systems interacting 
with high-dimensional environments will be increasingly 
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harder to distinguish. We can conclude that in such cases 
the effect of initial correlations has vanishing influence on 
the open-system evolution. 



IV. AVERAGE OPEN-SYSTEM TIME 
EVOLUTION: DEPOLARIZING CHANNEL 

In this section we consider the average time evolu- 
tion of an open quantum system. Making use of the 
decomposition of the unitary time evolution operator 
U t = W cxp{— iDt}W^ , we find for the average reduced 
system state 

( Ps (t)) = Ti E {(Wexp{-iDt}W^p(0)Wexp{iDt}W^)} 



{Pit)) 



(29) 



where p(0) denotes the initial state of system and envi- 
ronment. The averaged expression contains four unitary 
matrices and can therefore be obtained with the aid of 
the general result (A8), with Xi 
X 2 = p(0). This yields 



Xl = exp{— iDt} and 



(p(t)) = (W exp{-iDt}Wi p{0)W exp{iDt}W*) 



(f-d 2 \f(t)\ 2 
d 2 - 1 



I jd 



rf 2 l/WI 2 -i 

d 2 - 1 



p(0). (30) 



In quantum information theory this map p(0) — > (p{t)} 
is known as a depolarizing channel [16]. For the reduced 
dynamics we obtain the same map, 



(Ps(t)) = 



d 2 -d 2 \f(t)\< 
d 2 - 1 



I/d s 



d 2 \f(t)\ 2 
d 2 



1 



1 



Ps(0), 



(31) 



with the initial reduced state /?s(0) = Tr^p(0). Note 
that in contrast to the usual experience with open quan- 
tum systems, this map is always well-defined and com- 
pletely positive, depending only on the initial reduced 
state ps(0) regardless of the initial correlations which 
may be present in p(0) [22]. Further analysis may be 
carried out on the basis of this general result [12]. 



V. GENERIC DYNAMICS OF 
SYSTEM-ENVIRONMENT PURE STATE 
ENTANGLEMENT 

As another application of the results of Eqs. (7) and 
(18), we consider the special case of M — p, where p is 
an arbitrary fixed state of the total system. The term 
on the left-hand side of these equations then yields the 
unitary average value of the purity of the reduced system 
states. By Eq. (7) we obtain the time-independent value 
for a uniform average, to be discussed in the next section, 
while the result of Eq. (18) depends on time and on the 
eigenvalue distribution and will be discussed in Sec. VB. 




FIG. 3. (Color online) Average open-system purity within 
one standard deviation as a function of d_g for different values 
of d$- The asymptotic value is given by l/d$- Absolute and 
relative error decrease for increasing dimensions. 



A. The uniform average 

In the case of a uniform average over the unitary 
group we get the following connection between the av- 
erage value and the initial purity: 



(Hps)) 



(d 2 d E - d E )V(p) 



44 



i 



(32) 



It is remarkable that the right-hand side of Eq. (32) 
only depends on the purity P(p) of the initial state of 
the closed system and the dimensions of system and en- 
vironment. In the following we consider a pure state 
p = |*)(*| of the closed system with V{p) = 1. Then, 
the entanglement of the total state can immediately be 
read off the subsystem purity. A maximally mixed state 
in the subsystem corresponds to a maximally entangled 
state in the closed system, while a pure reduced state cor- 
responds to a separable total state [26, 27]. The variance, 
Eq. (12), can be simplified to 



Var[P(p s )] 



2(d|-l) 



(d s d E + l) 2 (d s d E + 2)(d s d E + 3) 
For large values of d E , the average purity behaves as 



(33) 



(rips)) 



i 



i - 



i 

4 



i 

d E 



(34) 



The asymptotic value 1/ds for d E — > oo corresponds to 
a maximally mixed state and, hence, to maximal amount 
of entanglement in the composite system. As it was 
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FIG. 4. (Color online) Generic subsystem purity evolution in 
a regular system for a pure initial state of the total system. 
All average values refer to Poissonian statistics. The generic 
evolution is plotted for two different systems, each for the 
examples of an initially pure and a maximally mixed state of 
the reduced system. After a few oscillations the asymptotic 
value is reached. The dependence of this value on the initial 
reduced purity vanishes for larger dimensions. 



foreshadowed by the general considerations displayed in 
Fig. 1 we find that the variance vanishes in the limit 
of large environmental dimensions, Var["P(ps)] — > for 
d,E —> oo. The same holds true for the relative fluctua- 
tions which for large cLe behave as 



y/Var[P(p S )] 

(Hps)) 



'2(4 



i) i 



4 



(35) 



For high dimensions the pure state entanglement between 
system and environment after the time evolution with 
respect to an average unitary operator demonstrates a 
universal behavior. The asymptotic value depends only 
on the dimensions of system and environment: For large 
d,E the total system approaches the maximally entangled 
state while the fluctuations around the average value de- 
crease to zero. This confirms the conclusion that generic 
states of high-dimensional systems contain large amounts 
of entanglement [8, 28-32]. Fig. 3 shows the dependence 
of the average purity within one standard deviation on 
the environmental dimension for four different system 
sizes. 



B. The average for general random matrix 
ensembles 



By inserting M = p = |*)(*| into Eq. (18), we get a 
time-dependent average of the reduced state's purity, 

(V(ps(t))}=c 1 (t) + c 2 (t) + (c 3 (t)+c 4 (t))(V(p s (0))), 

(36) 

which additionally depends on the initial subsystem pu- 
rity (P(ps{0)))' This stems from the terms ||Tr,gp|| 2 and 
||Tr5p|| 2 , which are equal for a pure state. The aver- 
age corresponds to an average over the eigenvectors using 



the Haar measure while the eigenvalue distribution is im- 
plicitly included in the constants ci(t), . . . , 04(f). Since 
a pure state remains pure under unitary evolution, the 
above equation characterizes the entanglement between 
system and environment for any time t. Hence, Eq. (36) 
can be interpreted as the generic system-environment 
pure state entanglement evolution. 

After additionally averaging the coefficients with re- 
spect to a certain eigenvalue distribution, symbolically 
represented by the subscript EVD, the subsystem purity 
evolution is given by 



(P(ps(t))) 



EVD 



d s d E + 1 
(d s + d E ){(g(t)) E vD - 



2d s d E mf(t) 2 rm}) 



EVD) 



(d s d E - l)(d s d E 

2d s d E m.f(t) 2 r(2t)}) 



EVD 



l)(d s d E + 
- (s , (*))evd 



(dsd E -l)(dsd E + 3) 



3) 

(P(Ps(0))), 
(37) 



with g(t) = A\f(t)\ 2 - |/(2i)| 2 - 44|/(t)| 4 . The aver- 
age values of these coefficients can now be obtained by 
inserting the respective average values of the functions of 
f(t). The derivations of these are given in Appendix B. 

We start with an analysis of regular systems, described 
by Poissonian level statistics. Inserting the correspond- 
ing averages into Eq. (37) yields the generic pure state en- 
tanglement evolution in a regular complex open system. 
This function is plotted in Fig. 4 for bipartite systems of 
equal sizes. We see that the reduced purity reaches an 
asymptotic value after some transient time t. This value 
depends on the initial subsystem purity but the depen- 
dence weakens for high values of rig. The corresponding 
limit yields the same result as for the uniform average: 



lim lim (V(p s (t)))poi 



1 



(38) 



Next, we investigate the generic evolution of the same 
quantity for typical chaotic systems, represented by the 
Gaussian unitary ensemble. In order to simplify calcu- 
lations, we employ an approximation for large values of 
d which is introduced in the Appendix B 2. Then, the 
asymptotic value for large times t becomes independent 
of the initial subsystem purity, 



lim (P(p S (t))) GVE 

t— >oo 



ds 



d s d L 



1 



(39) 



An additional limit for large dE leads to the same result 
as the Poissonian and uniform averages, 



lim lim (V(ps{t))) G VE = 

t— >oo d E — >oo ag 



(40) 



The dependence on the initial subsystem purity is eluci- 
dated in Fig. 5 for a high-dimensional system. Since the 
Hilbert space dimension is large, all initial states evolve 
towards the same asymptotic value close to the maxi- 
mally mixed state for both averages — a dependence on 
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FIG. 5. (Color online) Dependence of the subsystem purity 
evolution on the initial value. The plot shows a system with 
dimensions ds = 32 and cLe = 128 and five different equidis- 
tant values of the initial purity between a pure and a maxi- 
mally mixed state. Inlays: Magnification of the oscillations, 
which for the GUE are one order of magnitude smaller than 
for the Poissonian distribution. 



FIG. 6. (Color online) A close up comparison of the average 
purity evolution is shown for ds = 4 and different environ- 
mental dimensions. The continuous lines correspond to the 
initially pure reduced state and the dashed lines represent the 
evolution of the initially maximally mixed reduced state. 



A. Closed system thermalization 



the initial state can only be observed for small times. 
The overall shape of the two average evolutions is very 
similar although the oscillations are less pronounced for 
the GUE. 

A comparison of the generic evolution between Poisso- 
nian eigenvalue statistics and the GUE is shown in Fig. 6. 
The universality for large Hilbert space dimensions can 
be observed graphically: The gap between the asymp- 
totic values for t — > oo of two different initial states van- 
ishes as the environmental dimension increases. 



VI. EQUILIBRATION OF GENERIC 
QUANTUM SYSTEMS 

As a final application for the general unitary average 
values, we study the thermalization process of a com- 
plex quantum system - a problem which has attracted 
a significant amount of interest recently [11, 17, 33-37]. 
Before dealing with the evolution of a generic open sys- 
tem we start out with the simpler case of closed-system 
dynamics. Observation of the distance between the state 
at time t and the thermal equilibrium state allows to 
draw conclusions about the thermalization process. In 
this chapter we will derive explicit equations for the av- 
erage evolution of this quantity depending on system pa- 
rameters and the initial conditions. These can be used to 
infer whether or not a system will thermalize completely, 
i.e., the distance gets arbitrarily small for large times. 
In this case also the average time needed to reach the 
thermal equilibrium state up to a given proximity can be 
obtained. 



The equilibration process can be quantified by the dis- 
tance of the state p(t) = Utp(0)uj at time t to the 
thermal equilibrium state p G = e~@ H /Z, where Z = 
Tr{e~^ ff } and (3 = 1/kT with the Boltzmann constant 
k and the temperature T . As in previous sections, we 
rewrite the time evolution operator as Ut = W e~ tDt W^ . 
In this decomposition the thermal equilibrium state can 
be expressed as p G = W(e^ l3D ' /Z)W^ . For the Hilbert- 
Schmidt distance we get 

\\PG-P(t)\\ 2 

= \\W(e- 0D /Z)W^ - We~ tDt W^ p(0)We lDt W^\\ 2 
= \\W{e-P D /Z)W^f + \\We~ lDt W^p(0)We lDt Wm 2 

- 2Tr {W(e- 0D /Z)W^We~ lDt W^ p(0)We iDt W^} . 

(41) 

Integration over W employing the Haar measure yields 

(\\p G - P (t)\\ 2 )=V(p G ) + V(p(0)) 

- 2Tr{e-^ /Z (W^ p(0)W)} 

v v ' 

I/d 

= V(p G )+V( P (0))-2/d, (42) 

where we have used the invariance property of the trace 
under cyclic permutations to eliminate some of the uni- 
tary operators. Note that the result does not depend 
on time any more, even though the time dependence has 
been preserved in the unitary averaging process. It can 
immediately be confirmed that the value zero is reached 
for the case where both states p(0) and pa are maxi- 
mally mixed. This corresponds to a thermal state of in- 
finite temperature and minimal purity of 1/d. All other 
cases yield nonzero values since the purity is bounded by 
1/d < 'P(p) < 1. The purity of the Gibbs state can be 
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FIG. 7. (Color online) The dependence of the average thermal 
purity on the inverse temperature j3 is plotted for Poissonian 
level spacings and the GUE with ds = cLe = 2. Increasing 
temperature decreases the average purity as the Gibbs state 
gets closer to the maximal mixture. At low temperatures, 
regular systems have generically a higher thermal purity than 
chaotic systems. 



expressed in terms of the eigenvalues of H: 

V(pa) = Tr e - 2 ^/(Tr e -^) 2 = /(-2i/3)//(-i/ 



(43) 



with f{t) as in Eq. (19). Due to the negative exponents, 
the average value of V(pg) for different eigenvalue dis- 
tributions cannot be tackled with the same straightfor- 
ward approach which was used for the previously occur- 
ing functions of f(t). For this term we have to explicitly 
carry out the integration 



(P(pg)) = dEi... dE d P(E u ...,E, 



Eti ^ 
(44) 



where P(E\, . . . , Ed) denotes the probability for the set 
of energy levels {Ei, . . . , Ed}, see e.g. [15]. While for the 
Poissonian level spacing the energy levels are uniformly 
distributed, we get for the GUE [15] 



Pgue{Ei 



, Eh. 



Hi* 



Ej) 2 exp 



E^ 2 

i=l 



(45) 



with a normalization constant c. For both classes of 
systems this integral can be evaluated numerically for 
low dimensions. Furthermore, the absence of correlations 
in the regular spectrum allows to obtain average values 
of Eq. (43) even for higher dimensions by Monte Carlo 
sampling of equally distributed eigenenergies within the 




spectral span, which is determined by the normalization 
of the interaction strength. The dependence of (V(pg)) 
on j3 is shown for 2 x 2-systems in Fig. 7 for both distri- 
butions and Fig. 8 shows how the average thermal purity 
of regul ar systems depends on d E . 

T n I r I I 



10 20 30 40 50 60 70 80 90 100 

dE 

FIG. 8. (Color online) Average purity of the Gibbs state for 
regular systems with fixed values of ds = 2 and f3 — 10 as 
a function of dE- The data was obtained by Monte Carlo 
sampling of energy values and each point corresponds to an 
average over 10, 000 samples. 



B. Open system thermalization 

Next we want to refine our approach to the thermaliza- 
tion of quantum systems in contact with an environment 
[11, 37]. To this end we deal with the term 

\\Tr E {p G - P(t)}\\ 2 
= \\Tr E {W(e~ 13 D /Z)Wm\ 2 

+ \\Tr E {We- tDt W^ p(Q)We lDt W^}f 
- 2Tr {Tr E {W(e- 0D /Z)W^} 

x Tr E {We- lDt W^ p^We^W^}} . (46) 

After taking the unitary average of this quantity, we can 
identify the first two terms with already known expres- 
sions. Inserting M — e^ 130 jZ into Eq. (7) yields 

Tr E {W(e-? D /Z)Wl}\\ 2 ) = C\V(p G ) + C 2 , 

and from Eq. (18) we obtain for M = p(0), 

Tv E {We- lDt W^p{Q)We lDt W ] }\ r 



(47) 



= ci(t)P(/>(0)) + C2(t) + c 3 (t)V(ps(0)) + C 4 (i)P(pE(0)). 

(48) 

Averaging the third term requires knowledge of the sixth 
moment function of the unitary group: 
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2Tr {Tr E {W(e-P D j Z)W^x E {W e~ lDt W^ p(0)We iDt W^}} 
= 2^2 (<PiXj\ W(e-? D /Z)Wi \<p kXj ) {^kXi\ We' wt W^p(0)We tDt W^ \<p iXl ) 

ijkl 

= 2Y d ( ViXj \W{e-f >D /Z)Wi{I®\xi) (xi\)We- iDt Wip(0)We iDt Wi \<p iX i) ■ (49) 

ljl (...>^EC6)(e-^/ z ,Zig,| Xj >( Xi |,e-' D ',p(0), e * D ') 



The derivation of the general function E^ 6 ^ is given 
in the appendix A3. Here, we are particularly 
interested in the special case of E^(e~^ D /Z, I ® 
\Xj) {xi\ > e~ lDt , p(0), e lDt ), implying the choice of opera- 
tors 



X 1 = e-P D /Z,X 2 =I®\ Xj ){xi\, 



X* = XI 



-iDt 



X 4 = p(0), 



(50) 



which yields the simple result 



ijl 



x E< 6 ) (e-P D /Z, I ® \xi) (Xi\ , e- iDt ,p(0), e lDt ) \ ViX i) 
= Vd s . (51) 

Inserting this into Eq. (49), we finally obtain 



\Tt e { Pg - p(t)}\\ z ) = CxP{p G ) + C 2 + ~ Cl {t)V{ P m 

2 



C2(t) + c 3 (t)T( Ps (0)) + c^t)P(p E (0)) 



(52) 



The coefficients C\,C 2 and c\{t), . . . , 64(f) are well- 
known and the latter are also functions of f(i). Aver- 
aging this result for certain eigenvalue distributions of 
the Hamiltonian again enables us to study the influence 
of level statistics and quantum chaos [34]. We find re- 
markably simple expressions for the asymptotic values of 
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FIG. 9. (Color online) Generic equilibration in the limit of 
infinite-dimensional environments for regular and chaotic sys- 
tems. 



large d E : 



lim (\\Tv E {p G - p(t)}\\ > =c 



and 



lim (\\TT E {p G - p(t)}\\- ) =c 

ds-J-oo \ I GUE 



cos(i) sin(i) 



J x (2t) 



1 4 



t 



(53) 



(54) 



with Co = V(ps(0)) — 1/ds- We see that both functions 
describe an oscillating, algebraic decay to zero. Thus, 
for both eigenvalue distributions we find on average an 
asymptotic convergence to the thermal equilibrium state, 
see Fig. 9. While in the case of Poissonian level statistics 
the average squared Hilbert-Schmidt distance decays as 
<~ 4 , the same quantity behaves as t~ 6 for the Gaussian 
unitary ensemble. 



VII. CONCLUSIONS 

Employing general group theoretic results of Collins 
and Sniady [13] we have derived the explicit form of a 
class of operator-valued averages over the Haar measure. 
These were used to provide a tool-box for a variety of 
physically motivated applications. The thereby obtained 
expressions can be interpreted as averages over all pos- 
sible time evolutions. Motivated by the foundations of 
statistical mechanics, we assume that the dynamics of 
generic complex quantum systems is well represented by 
the average value. This is endorsed explicitly for high- 
dimensional systems by our results, showing a vanishing 
variance in the limit of d E —¥ 00 which is also in agree- 
ment with concentration of measure arguments. Due to a 
generalization of the uniform result for general eigenvalue 
distributions of the total system, detailed analysis of the 
effect of level statistics can be performed using methods 
of random matrix theory. 

To exemplify the general applicability and flexibility 
of our results we have investigated a series of topics as- 
sociated with open quantum systems and related fields, 
which are of current interest. We obtain particularly sim- 
ple results for the Hilbert-Schmidt distance of two quan- 
tum states in the open system, which was applied in the 
context of system-environment correlations. Our results 
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suggest an observable signature of the total state correla- 
tions on a generic open-system evolution, which vanishes 
in the limit of large cIe- 

Furthermore, we have investigated the generic evolu- 
tion of system-environment entanglement in a pure total 
state. Using exact analytical expressions, we have found 
a universal evolution towards the maximally entangled 
state with vanishing variance in the limit of large envi- 
ronmental dimensions. These results were confirmed by 
an analysis of the generic time evolution for regular and 
chaotic systems. 

Finally, we have studied the thermalization process 
of complex quantum systems. The average distance of 
the time-dependent reduced density operator to the re- 
duced thermal equilibrium state was given as a func- 
tion of time depending on system parameters and ini- 
tial conditions. While generic closed systems never com- 
pletely evolve from a non-equilibrium condition into equi- 
librium, high-dimensional environments lead to thermal- 
ization of generic open systems: As we have shown the 
average squared Hilbert-Schmidt distance exhibits an al- 
gebraically decaying oscillation, where however regular 
and chaotic systems decay with different powers of time. 

Using the results presented in this paper, expectation 
values for dynamical expressions of complex open sys- 
tems can be investigated analytically. This may be es- 
pecially useful if instead of the exemplary analysis of a 
particular model system, the expected performance of a 
method in generic systems is of interest. 
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Appendix A: Moments of the unitary group 

In this section we summarize our findings on the uni- 
tary averages based on the general results of Collins and 
Sniady [13]. For additional details and terminology we 
refer the reader to the standard literature on group the- 
ory, see, e.g., Refs. [38, 39]. Reference [13] provides a 
method to obtain exact expressions for average values of 
the form 

E<" +1 \Xi,...,X q )= [ dfi{U)UX 1 U^X 2 U ...UX q U\ 



(Al) 

where X\ , . . . , X q are arbitrary operators on a d- 
dimensional Hilbert space "H and dy, denotes the Haar 
measure on the unitary group U{d). The expression in- 
volves an integration over n = q + 1 unitary matrices 
and is therefore called an nth moment function of the 
unitary group. While for higher-order terms we have to 
make use of the above result, one can obtain moments 
up to the order of q = 4 also by elementary arguments 



using only the invariance property of the Haar measure, 
dfi(U) =dfi(UV) = d(j,{VU),V U,V e U{d). The second 
moment function E^(X) = J d(i(U)UXW is a trivial 
consequence of this invariance property. In fact, E( 2 )(AT) 
commutes with every unitary matrix V since 

VE<V{X)V* = [ du(UC/)UC/A(UC/) t 



= J dp,{W)WXW ] 



= E {2 \X). 

Using the irreducibility of the group of unitary matrices 
and Schur's Lemma we conclude that E< 2 ) (X) is propor- 
tional to the d-dimensional identity matrix: E( 2 )(X) = 
cl . The constant factor c can be obtained by taking 
the trace of both expressions: TrE( 2 )(A) = cd = TrA, 
and we obtain the general form of the second moment 
function as E' 2 )(X) = TrX • I/d. An elementary deriva- 
tion of the fourth moment function has been provided in 
Ref. [101. 



1. The general expression for the nth moment 

The elements of the matrix denoted by a capital let- 
ter will be denoted by the corresponding lowercase let- 
ter, leading to the correspondence of matrix elements xfj, 
Uij and matrices X a , U respectively. We rewrite the ex- 
pression for the nth moment function E^ 9+1 - 1 (X\, . . . , X q ) 
with n = q + 1 in terms of the elements: 



El 1 # q * 



'fc 2 „_i ki 



k 1 ,...,k 2 q 

x ( m kl . . . u k2q _ 2k2q _ 1 u ksk2 . . . Uj k2q 



(A2) 



u jk 2q > 18 



The expression (u ikl . . . u k2q _ 2k2q _ 1 u* k3k2 
called the nth moment of the unitary group, involving 
the average of a product of n elements of a unitary ma- 
trix. Since all odd moments are equal to zero [13], we 
write n = 2m. The result for the nth moment is a con- 
sequence of the Schur-Weyl duality [13]: 



ESi.i' ...Si i> 5i.ii ...Si ii Wsfxa -1 ) , 

<J,reS m 

(A3) 

where 



Ahr 



XA(e) 2 XA(g-) 



(A4) 
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1 


TABLE II. Number 


of constituents of the conjugacy 



TABLE I. Character table for the irreducible representa- classes of S3, 
tions of S3 labeled by partitions of 3. 



denotes the Weingarten function, and the Schur func- 
tion evaluated on d copies of x is written as s\.d{x) = 
s\ y d(x, . . . , x). The sum in Eq. (A4) is carried out over 
all partitions A of to, denoted by A h to. The charac- 
ters of the irreducible representations of S m labeled by 
partitions A and evaluated on the element a are repre- 
sented by xa (c) and e denotes the identity. The value of 
the Schur function for x = 1 can be determined with the 
help of the following lemma [13]: For any A h m one has 



«M(1) 



ml * — ^ 



(A5) 



which yields: 



EW(X U X 2 ,X 3 



dTrX 3 X 1 -TrAiTrA 3 



d{d 2 - 1) 
dTrX 1 TrX 3 -TrX 3 X 1 
d(<P - 1) 



(TrX 2 )J 



Ao 



(A8) 



An alternative derivation can be found in Ref. [10]. 



where c(r) denotes the number of cycles of r and Z T the 
number of elements for the conjugacy class associated 
with t. With this result it is possible to obtain the ex- 
act expression for the nth moment of the unitary group 
given that the character table of S m with respect to the 
irreducible characters labeled by partitions of to and the 
number of constituents of the corresponding conjugacy 
classes are known. 



2. The fourth moment 

The fourth moment contains an integral over four uni- 
tary matrices, hence, according to Eq. (A3), we have 
to = 2. The conjugacy classes and characters of the 
symmetric group S 2 are easily obtained since the group 
£2 on ly consists of two elements: The identity opera- 
tion, corresponding to the partition (1, 1), and the swap 
operation with the partition (2), each one forming its 
own conjugacy class. The two associated representations 
are the trivial representation and the sign representation. 
With Eqs. (A4) and (A5), the two Weingarten functions 
are 



Wg((l,l)) = ^- I , 
Wg((2)) 



d(d? - 1) 



(A6) 



With this we obtain the following average value with the 
aid of (A3), 



(e^(x 1 ,x 2 ,x 3 )) = 



T 1 T 2 
Xi,iX r 



kl x mn x op { u ik u no u m l u jp) , 

(A7) 



3. The sixth moment 

We derive the sixth moment function (to = 3) with the 
aid of the character table of S 3 with respect to the irre- 
ducible characters labeled by partitions of 3, see Tab. I, 
and the number of elements of all the three conjugacy 
classes, see Tab. II. The Weingarten functions for the 
three conjugacy classes of S3 are 



Wg((l,l,l)) = 
Wg((2,l)) = 
Wg((3)) 



-2 + d 2 



d(4 - hd 2 + d A ) ' 
1 

= -4 + 5d 2 - d 4 ' 
2 

" Ad - 5d 3 + d 5 ' 



(A9) 



The general sixth moment function is given in matrix 
elements as 



(E^(X 1 ,X 2 ,X 3 ,X i ,X 5 )) lJ 

= ({UX 1 U^X 2 UX 3 U^X i UX 5 U^)) lJ 

I \ "* 1*2 3*4 5*\ 

= ( y ^ UikX/ l .iU m iX mn 'U no X pU q pX qr U rs X st Uj t ) 

\k,...,t I 

= x kl x mn x\ p x qr x\ t (uikU n . u rs u m iU qp Uj t ') (A10) 
fc,...,t 



The final result can be obtained with the aid of Eq. (A3) 
and contains 36 terms, which we do not state explicitly. 
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TABLE III. Character table for the irreducible represen- 
tations of 54 labeled by partitions of 4. 



4. The eighth moment 

The eighth moment function (m = 4), expressed 
element-wise reads 



E^(X u ...,X r ) 



y t X kl X mn X op X qr X st X uv X wx^(a,T)^S ( Ta ) J 

(All) 



k,. . . ,x 



where 

d(cr.T) — $ia(rn)8na(q)dra(u)Ka(j)$kT(l)$oT(p)dsT(t)$wT(x)- 

(A12) 

The Weingarten functions for the five conjugacy classes 
of 1S4 can be obtained using Tabs. Ill and IV: 



Wg((4)) 
Wg((3,l)) 
Wg((2,2)) 
Wg((2,l,l)) 

Wg((l, 1,1,1)) 



ad ' 

-3 + 2d 2 
~a~d 2 : 
6 + d 2 



ad 2 ' 



-1 



9d - 10d 3 + d 5 ' 
6 - 8d 2 + d 4 



ad 2 



(A13) 



with a = -36 + 4M 2 - 14d 4 + d 6 . The sum of Eq. (All), 
generally containing 4! 2 = 576 terms, reduces to a sim- 
pler expression if the operators X\ , . . . , X7 are explicitly 
chosen for physical applications. 



Appendix B: Averages of functions of f(t) 

1. Regular systems 

Regular systems are characterized by Poissonian level 
spacings statistics [14, 15]. The energy levels are ran- 



domly distributed and show no correlations. In this ap- 
pendix we derive the average values of functions of the 
Fourier transform of the level density which continuously 
appear throughout the manuscript. We start with 



(i/wi 9 >=i(E 



-i{Ei-Ej)t 



d 2 

\ hi 
1 1 
d + tf 

id 



j dE x J dE 2 e-^- E ^ 



(Bl) 



R 2 (E U E 2 ) 



Since there are no correlations in the energy levels for 
Poissonian statistics, the 2-point correlation function is 
given by R 2 (E U E 2 ) = Ri(E 1 )R 1 (E 2 ). The level den- 
sity R\{E) is flat and homogeneous. Fixing the size 
of the matrix elements of the Hamiltonian by means of 
(\H tJ \ 2 ) = 1/d yields Rr(E) = ^{E + 2)6(2 - E) [9]. 
The resulting average of |/(t)| 2 is given by 



(l/(*)| 5 



'Poi — 



d- 1 

d 



sin(2i) 
2t 



(B2) 



and further average values, in particular those of the 
functions 3?{/(t) 2 /* (2*)} and \f(t)\ 4 can be derived by 
analogous methods. We obtain 
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<K{/(t)V*(2t)}) 



and 



(i/wi 4 ; 



Poi — 



f/ :! 



1 d- 1 

d 2 + ~dP~ 



3d -2 4(d-l) 



d 3 



sin(4t) 



4f 



,d-l 
J d 2 



sin(2i) 
2i 



sin(2i) 
2i 



2(d- l)(d-2) 
d3 



cos(2i) 



(d-2)(d- 1) 
d 2 



cos(2i) 



sin(2t) 
2i 



(B3) 



sin(2t) 



1 3 



2/ 



(d-l)(d-2)(d-3) 
d3 



sin(2i) 



1 4 



2f 



(B4) 



2. Chaotic systems: Gaussian unitary ensemble 

In this section we average the same functions of f(t) 
but now with respect to the GUE, representing chaotic 
systems with no time-reversal symmetry. We can build 
on the general results of the previous section using the 
n-point correlation functions for the GUE [15], 



R n (E u ...,E n ) =det 



where 



d~l 



.fe=0 



i,j=l ,n 

(B5) 



<p k {x) = (2 fc fc!y2Vd)- 1 /2 exp (-~*^ H k , 



(B6) 



and 



H k (x) = exp(x 2 ) (-—) exp(- 2 ; 2 ) (B7) 

denotes the Hermite polynomial of order k. The aver- 
age value of |/(t)| 2 is then obtained from Eq. (Bl). The 
average values of M{f(t) 2 f*(2t)} and |/(i)| 4 are deter- 
mined analogously. The integrations for the latter can 
be carried out with standard mathematical software but 
become fairly intricate for larger values of d. In fact, so 
as to facilitate computational analysis, it is possible to re- 
place (|/(i)| 4 ) ~ (|/(£)| 2 ) 2 - Correlations of fourth order 



are approximated by second order terms and the expres- 
sion can be considered as almost exact for our purposes. 
A more drastic approximation which allows analytical 
expressions for large values of d can be performed by ne- 
glecting all correlations and interchanging averages and 
powers of the function f(t) [9]. For the GUE this is valid 
in the leading order of 1/d [9]. Introducing h(t) = (f{t)) 
the three functions of interest are approximated by [9] 

(l/WI 2 ) * \Kt)\ 2 , 

(n{f(tfr(2t)})^di{h(t) 2 h*(2t)}, 

(|/(t)| 4 ) « \h(t)\\ (B8) 
The function h(t) is given by 



fc (*> = \ (E e ~ lElt ) = 2 1 dEe ^ Et V {6{E ~ Ej)) 



Ri(E) 



(B9) 



For the GUE this yields [9] 



d-l 



h G Mt) = \ J dEe-^J2^ k (E) Ji(2t)/t, 

(BIO) 



fe=0 



where J\ (t) denotes the Bessel function of the first kind. 
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